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ABSTRACT 

We combine the exceptional depth of the Hubble Ultra Deep Field (HUDF) images and the deep 
GRism ACS Program for Extragalactic Science (GRAPES) grism spectroscopy to explore the stellar 
populations of 34 bulges belonging to late-type galaxies at 0.8<z<1.3. The sample is selected based 
on the presence of a noticeable 4000 A break in their GRAPES spectra, and by visual inspection of 
the HUDF images. The HUDF images are used to measure bulge color and Sersic index. The narrow 
extraction of the GRAPES data around the galaxy center enables us to study the spectrum of the 
bulges in these late-type galaxies, minimizing the contamination from the disk of the galaxy. We use 
the low resolution (i?~50) spectral energy distribution (SED) around the 4000 A break to estimate 
redshifts and stellar ages. The SEDs are compared with models of galactic chemical evolution to 
determine the stellar mass, and to characterize the age distribution. We find that, (1) the average 
age of late-type bulges in our sample is ~1.3 Gyr with stellar masses in the range lO^'^-lO^" M©. (2) 
Late-type bulges are younger than early-type galaxies at similar redshifts and lack a trend of age with 
respect to redshift, suggesting a more extended period of star formation. (3) Bulges and inner disks 
in these late-type galaxies show similar stellar populations, and (4) late-type bulges are better fitted 
by exponential surface brightness profiles. The overall picture emerging from the GRAPES data is 
that, in late- type galaxies at z~\, bulges form through secular evolution and disks via an inside-out 
process. 

Subject headings: galaxies: spiral — galaxies: bulges — galaxies: formation — galaxies: evolution 



1. INTRODUCTION 

There are currently two alternative scenarios to ex- 
plain bulge formation in galaxies. First, semi-analytic 
models have traditionally proposed early formation from 
mergers, gen erating a scaled-down v ersion of an elliptical 
galaxy (e.g.. iKauffmann et al.lll993l ). Second, dynamical 
instabilities can co ntribute to the formation of a bulge in 
a primordial disk (jKormendv fc Kennicuttil2004f ). These 
instabilities can be triggered either i nternally, or by the 
accre tion of small satellite galaxies (jHernguist fc MihosI 
fT995h . and may result in later stages of star formation. 
Hence, the stellar populations in galaxy bulges provide 
valuable constraints to distinguish between these two sce- 
narios. 

The ability of the Hubble Space Telescope (HST) to 
resolve distant galaxies enabled t he study of bulges in 
galaxies out to redshift z :^ 1 (iBouwens et al.l 1l999l: 



Abraham etall Il999t lEUis et all 120011: iMenanteau et all 
2001; Ko o et alj 12003 : iMacArt hur et all I2008D . Simple 
p henomenologica l mod els — such as the one presented 
in IBouwens etal.1 (|1999| ) — have, in the past, tried to de- 
termine whether bulge formation happens before or after 
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the formation of the disk. The advantage of the lookback 
time probed out to z ~ 1 allows us to quantify the occur- 
rence of merging vs. secular formation of bulges. In the 
sample presented here we take advantage of the superb 
capabilities of the Advanced Camera for Surveys (ACS) 
to extract (slitless) low resolution spectra of bulges from 

faint galaxies at these redshifts.^ 

In their detailed review, iKormendv fc KennicuttI 
(|2004f ) discuss two distinct type of bulges, classical bulges 
— i.e., merger- built with Sersic index n > 2 — and 
pseudo bulges — built out of disk material having Sersic 
index n < 2. Early- type galaxies tend to have classi- 
cal bulges, while late-type galaxies are more likely to 
host a pseudo bulge. This scenario states that early- 
and late-type galaxies generally form their bulges in 
different ways. Many studies have been done at low 
and high redshifts to investigate properties of bulges 
and their formation histories. Studi e s on local galaxies 
(Ide Jond [l99l ICourteau et all [l99l iThomas fc DaviesI 
120061 : ICarollo etal.ll2007n have shown through colors and 
surface brightness profiles that later-type galaxies in the 
Hubble sequence have more bulges best-fit by an expo- 
nential profile ( disk-U ke) c ompared to an r^/^ profile. 
ICourteau et al.l ()1996D and Ide Jond (|1996D carried out 
bulge-disk decompositions for >80 galaxies, and found 
that 60-80% of late-type galaxies are b est-fit by the 
double exponential profiles. More recentlv. ICarollo et al.l 
((2007^ used HST ACS and the Near Infrared Camera 
and Multi Object Spectrometer (NICMOS) multi-band 
imaging to study the structure and the inner optical and 
near-infrared colors of local (z ~ 0) bulges in a sample 
of nine late-type spirals. Their analysis suggests that 
half of the late-type bulges in their sample must have 
developed after the formation of the disk, while for other 
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half, the bulk of stellar mass was produced at earlier 
epochs — as is found in early-type spheroids — and 
henc e must have developed b efore the formation of the 
disk. iThomas fc Daviei ()2006f ) analysed the central stel- 
lar populations of bulges in spiral galaxies with Hubble 
types Sa to Sbc by deriving luminosity- weighted ages and 
metallicities. They find that bulges are generally younger 
than early-type galaxies, because of their smaller masses. 
They suggest that bulges, like low-mass ellipticals, are 
rejuvenated, but not by secular evolution processes in- 
volving disk material. 

On the theoretical side, semi-analytical and iV-body 
simulations of galaxy formation have been mainly based 
on two basic assumptions. In the first scenario, all 
bulges result from the merging of disk galaxies (e.g., 
iKauffmann et al.l 119931 ). whereas the second one is 
based on an inside-o ut bulge formation scenario (e.g., 
Ivan den BoschI Il998f) . in which baryonic matter of a 
proto galaxy virializ e s and settles in an inside-out pro- 
cess. lAthanassoulal ()200 8^ argues that in order to ade- 
quately describe the formation and evolution of disk-like 
bulges, simulations should include gas, star formation 
and feedback. The author also states the importance 
of cosmologically-motivated initial conditions in the sim- 
ulations, since the properties of pre-existing disks may 
influence the properties of t he disk-like b u lges. When 
accounting for these effects, lAthanassoulal ()2008l ) simu- 
lated bulges that show properties similar to the observed 
disk-like bulges. 

The initial studie s of bulges at high redshift 
(lAbraham et al.lll999l : fEllis et al.ll200H iMenanteau et all 
l200lli " were done by measuring optical colors in galaxies 
to dis tinguish between bulge and disk colors. lEllis et aP 
()2001h analyzed the internal optical colors of early-type 
and spiral galaxies from the Hubble Deep Fields (HDF, 
IWiUiams ctal. 1996) for rcdshifts z < 0.6. They find 
that bulges are redder than the surrounding disks, but 
bluer than pure elliptica l s at th e same redshifts. In other 
work, IMenanteau et al.l ()2001l ) find strong variations in 
internal/central colors of more than 30% of the faint 
spheroidals in the HDF. They do not find such large 
variations in cluster galaxies, and hence estimate that at 
z ~ 1 , these strong color variations in field bulges are due 
to i nore recent episodes of star-formation Recent stud- 
ies (|Koo et al.ll200l iMacArthur et al.ll2008D have focused 
on bulge-disk decompositions to inve stigate the colors 
and radial profiles of bulges at z ~ 1. iKoo et all (|2005l ) 
present a candidate sample of luminous high-redshift 
(0.73 < z < 1.04) bulges (/sm < 23.1 mag) within the 
Groth Strip Survey, and find that majority of luminous 
bulges at z ~ 1 are very red. Their data favors an early 
bulge- formation scenario in which bu lges and field E-SO's 
form prior to disks. Mac Arthur et al.l (,2008 ) study bulges 
of spiral galaxies within the redshift range 0.1<z<1.2in 
the Great Observatories Origins Deep Survey (GOODS) 
fields, and find that bulges of similar mass follow similar 
evolutionary patterns. 

In this paper, we use the extraordinar y imaging depth 
of th e Hubble Ultra Deep Field (HUDF, iBeckwith et all 
l2006f ) and deep slitless grism spectroscopy using ACS 
from the GRism ACS Program for E xtragalactic Science 
(GRAPES) pro ject (PI: S. Malhotra; lElzkaLet al. 2004'; 
iMalhotra et al.L 20051 to explore the stellar ages of bulges 
in late- type galaxies at z ~ 1. The exceptional angular 



resolution and depth of the GRAPES/HUDF data com- 
bined with excellent grism sensitivity allows us to extract 
spectra of the most central regions of faint galaxies at 
z~l. 

This paper is organized as follows. The HST ACS ob- 
servational data, sample selection and photometric prop- 
erties of the selected late-type galaxies are presented in 
§[11 In § [3l we discuss morphological classification of 
these galaxies. The non-parametric CAS measurements 
are shown in §12111 while we use GALFIT in l3.2l to an- 
alyze galaxy morphologies in two dimensions. We study 
the bulges stellar populations via Spectral Energy Dis- 
tribution (SED) fitting of their GRAPES spectra in § H 
Our results are discussed in § together with the pos- 
sible biases due to the age-metallicity relation, and our 
conclusions are summarized in § |6l 

Throughout this paper we refer to the HST ACS 
F435W, F606W, F775W, and F850LP fihers as the B-, 
V-, i'-, and z'-bands, respectively. We adopt a Hubble 
constant Ho=70 km s~^ Mpc~^ and a fiat cosmology 
with fim=0.3 and ^1a=0.7. At redshift z ~ 1, this cos- 
mology yields a scale of 1"— 8.0 kpc. The lookback time 
is 7.7 Gyr for a universe that is 13.5 Gyr old. Ma gnitudes 
are given m the AB system ()Oke fc Gunnlll983[ ). 

2. OBSERVATIONAL DATA AND SAMPLE PROPERTIES 

2.1. The UST/ACS Data 

The HUDF is a 400 orbit survey of a 3.4' x 3.4' field 
carried out with the ACS in the B, V, i' and z' filters 
(|Beckwith et al.ll2"006f ). We have carried out deep shtless 
spectroscopy of this field with the ACS grism as a part 
of the GRAPES project, which was awarded 40 HST or- 
bits during Cycle 12 (ID 9793; PI: S. Malhotra). The 
grism observations were taken at five different orienta- 
tions, in order to minimize the effects of contamination 
and overlapping from nearby objects. The details of the 
observations, data re duction and the fina l GRAPES cat- 
alog are described in lPirzkal et al.l (2004), who extracted 
the grism spectra for objects in the HUDF to a limiting 
magnitude of z^g ~27.5 mag. These spectra cover the 
wavelength range between 6000 to 9500 A with a resolv- 
ing power of i? ~ 50 — 100 ^ and are characte rized by a 
net significance 7V>5 fsee iPirzkal et al.ll2004D . We used 
the multi-band high-resolution HUDF images to study 
each object at z ~ 1 in the ACS band closest to the B- 
band rest-frame, in order to min imize any effects of th e 
morphological K-correction fe.g. JWindhorst etal.ll200"^ . 

2.2. Sample Selection and Properties 

The ACS grism sensitivity peaks at ~8000 A, and 
hence the GRAPES spectra are sensitive in identifying 
galaxies at z~ 1 through their 4000 A breaks, and galax- 
ies at z~5 — 6 through their Lyman breaks. We selected 
objects with high signal-to-noise in the one-dimensional 
(ID) GRAPES grism spectra (-1500 in the HUDF). 
All these spectra were extracted using narrow extrac- 
tion windows (5 pixels wide — the pixel-scale in the 
grism images is 0.05"/pixel) around the center of each 
galaxy. We compared the running average of fiux in 
10 data points between neighboring regions on each ID 

* Slitless spectroscopy produces a variable spectral resolution, 
depending on the siz e of the object . The d etails concerning this 
issue are described in lPasquali etafl l l2006al ) 
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Fig . 1 . — Distr ibution of amplitude of 4000 A break (D4000; 
IBaloeh et aI|[T999t) for all gal axies in our sample (grey histogram). 
D4000 was measured as a ratio of average continuum flux redward 
and blueward of the 4000 A break from ID GRAPES grism spec- 
tra. The average value of the D4000 for our samp le is ^X.'Si and for 
com parison, elliptical galaxies have D4000 > 1.6 f Kaufl^ma nn et ahl 
|2003; Padmanabhan et al. 2004). We have also plotted D4000 
(hashed histogram) for GRAPE S /HUDF early- type galaxies (0.6 < 
2<1.1). l lPasauah et al.ll2006bl') . 

spectrum. If the difference between two neighboring, av- 
erage flux values was greater than Sct, we considered that 
change in flux level as a 'break' in the spectrum. A large 
('^ISO) number of objects that show continuum breaks 
were selected using this technique. After visual selection 
procedure, many objects were cla ssified and studied as 
Lyman break gala xies at z > 4.5 (jMalhotra et al.l l2005t 



Hathi et all [20b8al) . or as late-ty pe stars (jPirzkal et alj 
2005D, or as ellipticals at z ~ 1 (|Pasauah et al.l l2006bD . 



During this visual classification, we also found that a 
number of late-type galaxies (mostly spirals) showed a 
prominent 4000 A break (at observed ~8000 A), which 
is the major spectral feature due to the presence of an 
old stellar population. This feature was in general ob- 
served in all grism spectra obtained for each object at 
different position angles and for which spectral contam- 
ination was negligible. We selected 34 late-type galaxies 
(median redshift 2:~ 1), to study the properties of their 
central stellar populations. We me asure D4000 — the 
amplitude of the 4000 A break (e.g.. iBalogh etal.lll999h 
— as the ratio of the average continuum flux redward and 
blueward of the 4000 A break from the ID grism spec- 
tra. Figure [1] shows the distribution of D4000 for these 
galaxies at z~l. The average D4000 (^1.3) for the se- 
lected late-type galaxies (grey histogram) is smaller than 
the D4000 observed for a typical e lliptical galaxy (> 1.6, 
iKauffmann et all 120031 : fPadmanabhan et all I2004f ). We 
also show D4000 (hashed histogram) for e lliptica l galax- 
ies (0.6 <z< 1.1) from iPasauah et all (|2006b[ ). Fig- 
ure [2] shows the color-composite images of 6 represen- 
tative late-type spiral galaxies at z ~ 1 in our sample, 
which clearly demonstrates that all have redder bulges 
in their centers. 

Table [T] shows the optical [BVi'z') magnitudes, 
the HUDF IDs, and coordinates for each selected 
galaxy, as obtained f rom the published HUDF catalog 
(jBeckwith et al.ll2006[ ). The last column of Table [1] gives 
the observed iy-z'~) colors for these galaxies, correspond- 
ing to roughly rest-frame {U-B) colors at z~l. Table[T] 




Fig. 2. — Color composite images of a representative sample of 
late- type spiral galaxies at 2 ~ 1 . Axes show size of the stamp in 
pixels (0.03"/pixel). Note that all have small central bulges that 
are (in general) redder than their disks. 

gives all available redshifts for these galaxies. The first 3 
columns in Table [T] show the published HUDF IDs and 
coordinates, while the fourth column shows the photo- 
metric redshifts from the GRAPES sp ectro-photometric 
redshift catalog (for description, see iRvan et all l2007| ) 
or the GOODS-MUSIC (MUltiwavelength S outhern In- 
frared Catalog) catalog ( Grazian et al.l 1200^ . The fifth 



colum n of Table [T] gives th e spectroscopic redshi fts from 
VLT (|Grazian et al.l l2006t IVanzella et cdl l2008f ). when 
available. The last column of Table [T] gives the redshifts 
from GRAPES SED fitting as described in § H The 
latter are the redshifts used throughout this paper. No- 
tice that estimating the redshift of faint spectra without 
emission lines can be challenging even for deep surveys 
such as VVDS (Le Fevre ctal. 2005). Our slitless grism 
data have an optimal spectral resolution that maximises 
the S /N for populations with a prominent 4000 A break. 
Comparisons of redshift estimates of faint (iab —22-24 
mag) early-type galaxies at z ~ 1 from the Probing Evo- 
lution And Reionization Spectroscopically (PEARS) sur- 
vey show that ACS slitless grism data often fare better 
than publicly available spectroscopic redshifts (Ferreras 
et al. 2008, in preparation). 

2.3. Observed Color Profiles 

We used the Interactive Data Language (IDL^) pro- 
cedure APER^" to compute aperture photometry using 
several aperture radii. We chose our starting aperture 
radius to be 2.5 ACS pixels (^O'.'l), because the width 
of the narrow extraction window (see lPirzkal et al.|l2004l . 
for extraction details) of our GRAPES spectra is 5 pixels. 
We measured aperture magnitudes for all galaxies in our 
sample in all four ACS bands {BVi'z') available from the 
HUDF with aperture radii ranging from 2.5 pixels to 35 
pixels. Using these aperture magnitudes, we measured 
the {V-z') color profiles for our galaxies. Figure [3] shows 
these color profiles for 6 sample galaxies. This figure 
show that the inner disk in most of these galaxies is red, 
and is dominated by the older stellar population. Our 
galaxies show that for all apertures the {V-z') color is 
redder than 0.9, with redder colors at smaller aperture 
sizes (the central bulge region). 

Figure S] shows the observed {V-i') color as function 
of redshift for bulges and spheroids. The black dots cor- 
respond to our sample of HUDF/GRAPES bulges. The 
crosses are the GOODS/HDF-N spheroids presented by 



^ IDL Website |http: //www. ittvis. com/index. asp| 

IDL Astronomy User's Library 

|http: / /idlastro.gsfc.nasa.gov/homepage.html| 



Website 
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Fig. 3. — Observed (V-z') colors measured using aperture mag- 
nitudes in eight different size apertures. Here we show colors for 
six sample galaxies. These are the same six galaxies shown in Fig- 
ure |2] Note that with the exception of object ID 501, all these 
galaxies become bluer from the inside outwards. 

iMacArthur et al.' (2008), and the filled triangles are the 
GOOD S /CDF-S early- type galaxies from Ferrer as et al.l 
()2005f l. The lines represent the expected color evolution 
for a set of star forma tion histories using the p opula- 
tion synthesis models of iBruzual fc CharlotI (|2003D . The 
thick lines track the color evolution for a stellar popula- 
tion with solar metallicity and an exponentially decaying 
star-formation rate, starting at z^? = 5, with an exponen- 
tial decay timescale of 0.5 (solid), 1 (dashed) and 8 Gyr 
(dotted). The thin solid line shows the expected color 
evolution for a decay timescale of 0.5 Gyr and 1 /3 of solar 
metallicity. Most of the bulges in early-type galaxies in 
GOODS/CDF-S (|Ferreras et all [20051 : IMacArthur et all 
l2008f ) are consistent with short decay timescales, while 
the bulges in our late-type spiral sample agree better 
with a more extended star formation history. This figure 
shows that the bulges explored in this paper belong to a 
similar popula t ion as the "blue early-types" presented in 
iFerreras et al.l ()2005D . This blue population constitutes 
~20% of the total sample of early-type systems visu- 
ally selected in GOODS/CDF-S. Because of the excellent 
photometric depth of the HUDF, our sample ex tends the 
bulge redshift distribution of IMacArthur etal.l l|2008( ) to 
z~ 1.3. 

3. MORPHOLOGICAL PROPERTIES 

We performed a morphological analysis of the sam- 
ple galaxies in two steps: (a) a non-parametric analysis 
of the distribution of the galaxy light, using the mea- 
sures of asymmetry, concentration and dumpiness (or 
smoothness) to confirm our visual inspection; (b) a two 
dimen sional (2D) dec omposition performed with GAL- 
FIT (jPeng et al.l l2002l ) . in order to quantify the galaxy 
morphology and in particular to extract their Sersic in- 
dices. 

3.1. CAS Measurements 

We use the classical Concentration (C), Asymmetry 
(A) a nd clumpinesS (o r smoothness) — the CAS param- 
eters (jCons elice etaL| [2000: ConseUciHool) — to carry 
out the non-parametric approach to quantify morphol- 
ogy. We computed the C and A values for our galax- 
ies followijig_tlie de&^ and methods as discussed 
in IConselice et all ()2000[ ). The Concentration index cor- 
relates with the Sersic index and the Bulgc-to-Disk ra- 
tio: high C-values correspond to early-type morphology, 
while lower C-values are suggestive of a disk-dominated 
or later-type and irregular galaxies. Asymmetry can dis- 
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Fig. 4. — Color vs. redshift diagram of bulges and spheroids. The 
black dots correspond to our sample of HUDF/GRAPES bulges. 
The crosses are the GOODS/HDF-N spheroids presented by 
IMacArthur et al.l l|^08!), and the triangles are the GOODS/CDF-S 
early-type galaxies from Ferreras ctal. (2005). The lines represent 
the expected color evolution for a set of star-formation histories 
from the population synthesis models ofBruzual & Chariot (2003). 
The thick lines track the color evolution of a stellar population with 
solar metallicity and an exponentially decaying star-formation rate 
starting at zp = 5, with an exponential decay timescale of 0.5 
(solid), 1 (dashed) and 8 Gyr (dotted). The thin solid line shows 
the evolution with a timescale of 0.5 Gyr at 1/3 of solar metallicity. 

tinguish irregular galaxies or perturbed spirals from re- 
laxed systems, such as E/SO and normal spirals, dumpi- 
ness quantifies the degree of structure on small scale, 
and roughly correlates with the rate of star-formation. 
We derived the Concentration and Asymmetry indices 
using the images taken in the z'-band, which roughly 
corresponds to the rest-frame _B-band at z ~ 1. Our 
measurements of C and A are shown in Figure [5l to- 
gether with the mean loci for n earby early-, m i d- an d 
late- type galaxies as derived by iBershadv et al] (|2000D . 
A few perturbed galaxies show higher asymmetry val- 
ues. Figure [5] clearly shows that our sample of galaxies 
consists of mostly late-type galaxies, and confirms our 
visual morphological classification of late-type galaxies 
with bulges. 

3.2. 2D Galaxy Fitting using GALFIT 

GALFIT (jPeng et al.ll2002D is an automated algorithm 
to extract structural parameters from galaxy images by 
fitting/decomposing these with one or more analytic 2D 
functions. It offers different parametric models (the 
"Nuker" law, the generalized Sersic-de Vaucoulcurs pro- 
file, the exponential disk, and Gaussian or Moffat func- 
tions), and allows multi-component fitting, which is use- 
ful to measure Bulge-to-Disk (B/D) or Bulge-to- Total 
(B/T) hght ratios. 

3.2.1. Thumbnail Image Extraction 

The GALFIT disk-|-bulge decompositions were per- 
formed on thumbnail (or "postage stamp") images ex- 
tracted around the objects in our sample, rather than 
on the entire science image itself. Three thumbnail im- 
ages for each object were extracted from the original 
HUDF images. All thumbnail images are 201x201 pixels 
(^ 6'.'0x6'.'0) in size. The first thumbnail was extracted 
from the science image itself. The second thumbnail was 
extracted from the comprehensive segmentation image 
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Fig. 5. — Asymmetry and Concentration values for the selected 
late-type galaxies. The distinction between early-, mid- a nd late- 
type galaxies is from Bcrshady et al. (2000) and Conseli ce et al.l 
l|2005l ) . According to this classification scheme, our sample - se- 
lected by the presence of a prominent 4000 A break and removing 
the spheroidal galaxies — corresponds to late-type spirals. 

generated by ICoe et all ()2006f ). which was used as the 
"bad pixel map/mask" image used in GALFIT. GAL- 
FIT uses this "mask" image so that all non-zero val- 
ued pixels are ignored in the fit. Hence, the extracted 
segmentation stamps were modified, so that only pix- 
els belonging to the galaxy had zero value, while any 
pixels belonging to another object are set to a non- 
zero value. We tested this GALFIT decomposition with 
and without bad pixel maps for comparison, and ob- 
tained very similar fitting results. The third thumbnail 
was extracted from the drizzle-generated weight image 
iKoekemoer etal. 2002). These weight images were mod- 
ified for GALFIT (C. Peng, private communication) as 
follows. First, the science image was smoothed by a few 
pixels to get rid of some of the random pixel-to-pixel 
variations. Second, a variance image S, was calculated 
using S — (1/wht) + data/ exptime, where, wht is the 
drizzle-generated weight image, data is the science im- 
age in counts/sec and exptime is the total exposure time 
of the image. Finally, a sigma image is generated using 
(7 = ^/S. These modified weight images were used as 
the input "sigma" image (noise maps) in the GALFIT, 
which is necessary for proper error-propagation. 

3.2.2. Sky Background 

The drizzled HUDF images are sky subtracted and 
therefore, to understand the effects of sky-subtraction, 
we used the following procedure. A careful analysis of 
the HUDF sky-background and its corresponding uncer- 
tainties was performed by Hath ictal. ( 2008 b). We used 
the sky-background values from lHathi et al.l |2008bl ). and 
allowed GALFIT to either vary this sky-level during the 
fitting process, or keep it fixed. For comparison, we 
also used the sky-background measured from each in- 
dividual object stamp, and repeated the process. We 
also tested our GALFIT decomposition by comparing 
the results from GALFIT with and without the addition 
of the sky-background. We found very consistent and 
similar fitting parameters from all these tests. There- 
fore, we adopted th e sky-background levels measured by 
iHathi et all pOOSbl ) in our final GALFIT decomposition. 
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Fig. 6. — The distribution of Sersic indices in) for the galaxy 
retrieved from GALFIT single component Sersic fit to galaxies in 
V- and z'-band. The mean values of n is reported in top right 
corner. Using a two-sided K-S test on these distributions, we can- 
not reject the hypothesis that the V- and z'-band distributions are 
drawn from the same population. 

3.2.3. Using GALFIT 

GALFIT produces model images of galaxies based on 
initial input parameters. These images are convolved 
with the ACS Point Spread Function (PSF) image be- 
fore comparing with the actual galaxy image. Fitting 
proceeds iteratively until convergence is achieved, which 
normally occurs when the does not change by more 
than 5 parts in 10* for successive five iterations (see 
iPeng et al.ll200l for details). 

GALFIT requires initial gue sses for the fitting pa- 
rameters. Following |D ong & D e RobertisI (|2006f ) and 
ISimien k, de VaucouleursI (ll986) , we used the output pa- 
rameters from the pu blished HUDF SExtractor catalogs 
(|Beckwith et al.|[2006[ ) as input values for the magnitude, 
the half-light radius, the position angle, and the elliptic- 
ity of each object. The initial value for the Sersic index, 
n, was taken to be 1.5 (jCoe et al.ll2006[ ). Tests based on 
an adopted initial value of n=4 showed similar results. 

All GALFIT measurements were obtained from the V- 
and z'-band images (approximately rest-frame U- and 
B-band at z~ 1, respectively). First, we fitted only an 
one-component Sersic profile to our galaxies to improve 
the initial estimates for the Sersic index, the axis ratio 
and the position angle of each galaxy. The distribution 
of Sersic indices for our galaxies in the V- and z'-band is 
shown in Figure [S] Our measurements of the Sersic index 
in V- and z'-ban d arc comparable to i'-band values of the 
galaxies in lCoe et al. (2006). Figure [5] confirms that most 
of the galaxies in our sample have a Sersic index n < 2 in 
both V- and z'-bands, which implies that our galaxies are 
disk-dominated, as Figure [2] and Figure [3] clearly suggest. 

Next, we simultaneously fitted two components: a 
Sersic profile plus an exponential disk profile, to get bet- 
ter estimates of the bulge and disk magnitudes, respec- 
tively. For this simultaneous fit, we kept the coordinates 
of the galactic center (within ±1 pixel constraint), the 
axial ratio and the position angle fixed, while we al- 
lowed GALFIT to fit rribuige, fndisk, Re, Rs and the the 
Sersic index (n). For all galaxies in our sample, we ob- 
tained better fits for bulge Sersic indices of 71 < 1.0. We 
also tested our runs by fixing initial value of n=l for all 
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galaxies and found that GALFIT converged to similar 
solutions in the end. The bulge and disk models ob- 
tained from these best fits were then used to estimate 
the B/D ratio. We used an aperture of 5 pixels diameter 
to measure the B/D and B/T ratios, which used the same 
aperture size as for extracting the GRAPES SEDs. The 
majority of our galaxies show a B/D value <1 within this 
aperture. For larger apertures encompassing the galax- 
ies' total light, B/D appears to be << 1, in agreement 
with our galaxies being disk-dominated (i.e., late- type 
galaxies) . 

4. STELLAR POPULATION MODELS - BULGE 
PROPERTIES 

The GRAPES grism spectra were taken at five differ- 
ent position angles (PAs) to remove any contamination 
and overlap from nearby objects. We generated one fi- 
nal spectrum for each galaxy by combining all of the 
GRAPES spectra obtained at the 5 different PAs. The 
combination was performed as a simple averaging opera- 
tion, after resampling the spectra onto a common wave- 
length grid. Portions of spectra which were contaminated 
more than 25% (see iPirzkal etal..,2004 . for a description) 
were not used, unless absolutely necessary. The Poisson 
errors were propagated, and the standard deviation of 
the mean between the 5 individual PAs was computed. 
The larger of either the Poisson noise or the standard de- 
viation of the mean was used in the subsequent analysis. 
Our goal is to fit stellar population models to the age- 
sensitive 4000 A break observed in the GRAPES spectra 
of these galaxies. 

4.1. Star- Formation Histories (SFH) 

We fit our ACS grism spectra to a grid of models ob- 
tained by combining the simple stellar populations of 
iBruzual fc CharlotI ()2003[ ). A standard method is 
used. We explore a wide volume of parameter space 
in order to infer robust constraints on the possible ages 
and metallicities of the stellar populations in the central 
bulges of these galaxies. This comparison requires a care- 
ful process of degrading the synthetic SED (resolution 
R ~ 2000) to the (variable) resolution of the GRAPES 
spectra. Special care must be taken with respect to the 
change of the Line Spread Function (LSF) with wave- 
length, which results in both an effective degradation of 
the spectral resolution as a function of wavelength, and a 
different net spectral resolution with respect to the size of 
the galaxy. After exploring a range of values R = 30 — 80, 
we find that an effective resolution of i? = 50 is suitable 
for all the spectra in our sample. The ACS grism spectral 
resolution is not degenerate with respect to parameters 
describing the star-formation history, and mostly results 
in a global shift of the likelihood. 

In order to determine the redshift as accurately as pos- 
sible, we start with some guessing values obtained from 
three sources: a photometric redshift; a VLT spectro- 
scopic redshift — where available — and a redshift esti- 
mate taken from an automated method (as discussed in 
S I2T21) to search for a prominent 4000 A break in GRAPES 
data. A small set of templates at the GRAPES resolu- 
tion were used to determine the best redshift for each 
galaxy, using the guessing values described above as a 
starting point, and performing a simple cross-correlation 
for a range of redshifts until the best match is found. 



This method generates the redshifts used throughout this 
paper, shown in Tabled] as "SED Fit" (last column). 

In order to make a robust assessment of the ages and 
metallicities of the unresolved stellar populations, we use 
two different sets of models to describe the build up of the 
stellar component. The models depend on a reduced set 
of parameters, which can characterize a star-formation 
history in a generic way. 

Model #1 (EXP): We take a simple exponentially de- 
caying star formation rate, so that each star-formation 
history is well parametrized by a formation epoch, which 
can be described by a formation epoch (t(zi?)); a star- 
formation timescale (r^, = 0.1 ^ 4 Gyr); and a metal- 
licity {[m/H] — —1.5 -1-0.3), which is kept fixed 
at all times. The numbers in brackets give the range 
explored in the analysis of the model likelihood. The 
range in formation epochs is chosen from zp = 10 to 
t{zp) — 0.2 Gyr (this range depends on the observed 
redshift of the galaxy). 

Model #2 (CSP): We follow a consistent chemical 
enrichment code as described in Ferrer a s fc Silk (2000) . 
The model allows for gas infall and outflow. The metal- 
licity evolves accor ding to these para meters, using the 
stellar yields from iThielemann et all ([l9 96) for mas- 
sive stars (> 1071^0). and lvan den Hoek fc GroenewegenI 
(1997) for intermediate mass stars. The free parameters 
are the formation epoch (same range as the one cho- 
sen for Model #1), the timescale for the infall of gas 
(r/ = 0.1 ^ 1 Gyr), and the fraction of gas ejected in 
outflows (-Bout = — > 1). The star-formation efficiency 
is kept at a high valu e Ceff = 20 as expec ted for early- 
type populations (see iFerreras fc SiI3l2000f ). 

For each of the two sets of models we run a grid of 
SFHs, con volving simple stel l ar pop ulations from the 
models of IBruzual fc CharlotI (|2003l ). The grid spans 
64 X 64 X 64 SFHs (note that three free parameters are 
chosen in each set) over a wide range of values as shown 
above. Once the best fit is obtained within the grid, 
we run a number of models with random values of the 
parameters with an accept /reject criterion based on the 
likelihood - analogous to the Metropolis algorithm, e.g., 
ISahal (l2003h . The process ends when 10,000 models are 
accepted. The total number of accepted models deter- 
mine the median and confidence levels of the parameters. 
The distribution of reduced values has an average of 
X^=0.74 and RMS a{x^) = 0.35. The used through- 
out includes a mild Gaussian prior on the mctallicity with 
average [m/H]= -0.1 and RMS a{[m/H]) = 0.5. This 
prior allows for a relatively wide range of average metal- 
licities, and is compatible w ith the values com monly 
found in these systems (e.g.. ICarollo et all l2007t l. For 
the interested reader, we include an appendix describing 
in detail the effect of the application of this prior. We 
find no significant change in the estimates of the stellar 
age distribution with or without priors. 

Figure [7| shows the best fit models and the observed 
SEDs for 10 galaxies in our sample. The error bars rep- 
resent the observations, and the solid line corresponds 
to the best fits for the CSP model. The wavelength is 
shown in the observed frame, whereas the wavelength 
range chosen for all galaxies is 3800-5000 A in the rest- 
frame. This choice ensures a consistency in the com- 
parison of the stellar populations in our sample. The 
chosen range straddles the age-sensitive 4000 A break. 
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Fig. 7. — Spectral energy distributions of 10 galaxies in our 
sample (each one is labeled with the HUDF numbers and model 
redshifts). The error bars are the observed ACS/G800L data and 
the lines are the best fits according to the CSP models (see text 
for details). The distribution of reduced values for all galaxies 
in our sample has an averag e of x^=0.74 and RMS o-(x^) = 0.35. 
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Fig. 8. — The likelihood distribution with respect to average age 
for all 34 bulges. 

To illustrate the uncertainties of estimating stellar ages 
within the modelling, we show in Figure [5] the likelihood 
distribution with respect to average stellar age for each 
galaxy. 

Figure M shows the 4000 A break amplitude, D4000, 
as a function of stellar age. The distribution of bulge 
ages is overplotted as a histogram, which agrees very 
well with the observed range in D4000 - which is 1.2 
to 1.5, as shown in Figure [1] Different curves show 
three simple stellar population models with three dif- 
ferent metallicities. For metallicity around solar (thick- 
solid and dashed lines), representative of our bulges, the 
variation in D4000 with age is very similar. For solar 
metallicity with 'E{B-V)^0.2 dust reddening (thin-solid 
line) and for high metallicity (dotted line), the relation 
is somewhat steeper. Figure [9] shows that for the range 
in D4000 and metallicities of our sample, the effect of 
metallicity on D4000 is small, hence D4000 is a good age 
indicator for this sample. 

Figure (TU] shows the ages and metallicities of the best 
SFHs for each bulge. The average and RMS scatter for 
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Fig. 9.— The 4000 A break amplitude, D4000, as a function 
of stellar age. The distribution of bulge ages is overplotted as 
histogram, which agrees very well with the range in D4000 (see 
Figure [TJ. Different curves show three simple stellar populations 
model with three different metallicities, as labelled. The effect of 
dust on a solar metallicity model is shown as a thin solid line. The 
figure illustrates that for the range of ages found in our sample, 
the degeneracy caused by dust and metallicity will not result in a 
significant change of the average stellar ages. Furthermore, if the 
bulges were dusty or more metal rich, the resulting populations 
would become even younger, ruling out the possibility of bulges as 
old as early- type galaxies at the same redshift. 

age and metallicity are shown as dots and error bars, re- 
spectively. For the EXP models - which have zero spread 
in metallicity - the error bars in metallicity represent 
the uncertainty estimated from the likelihood. The solid 
lines in the lower panels correspond to the age of the Uni- 
verse as a function of redshift. The dashed lines show the 
age that a simple stellar population - a population with 
a single age - would have if formed at redshifts (from 
top to bottom) zp = {5, 3, 2}. The median value of the 
stellar ages of our bulges is 1.3 Gyr. The CSP models 
treat chemical enrichment in a more consistent way than 
the EXP models and should therefore better reflect the 
true populations. We use mass-weighted average ages 
for these stellar populations from the SFHs because they 
better reflect the formation process of bulges (or galaxies 
in general). A very small amount of young stars - some- 
thing that may not reflect the true formation process of 
the bulge - can have a large effect on luminosity- weighted 
ages. By using composite models such as EXP or CSP we 
minimise this contamination by using the mass-weighted 
ages instead. 

We would emphasize here that it is the average stel- 
lar age that can be reasonably constrained with the data. 
Therefore, FigurefTOldoes not imply that all stars in these 
bulges are ^--^1.3 Gyr but many may be older. To clarify, 
the formation epoch (characterized by a formation red- 
shift) is the age when star formation starts in the model. 
Comparing observations and models of unresolved stellar 
populations can only give us robust average ages (the first 
order moment of the age distribution) and, with higher 
uncertainty, we can also determine the "width" of the 
age distribution (the second order moment). We caution 
the reader that the actual parameters used in the mod- 
elling (especially formation redshift) constitute a way to 
characterise a generic set of star formation histories, but 
the uncertainties in these parameters will be larger than 
those in the average age, which we consider to be the 
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Fig. 10. — Ages and metallicities corresponding to the best fit ac- 
cording to a simple exponentially decaying model (EXP; [Left]) or 
a consistent chemical enrichment code (CSP; [Right]). The filled 
circles are the average values of age or metallicity and the error 
bars represent the RMS of the distribution. The solid lines in the 
bottom panel track the age of the Universe at a given redshift for a 
concordance cosmology. The dashed lines — from top to bottom - 
correspond to the ages that simple stellar populations would have 
if formed at redshifts zp = {5,3,2}. The grey triangles are the 
GRAPES/HUDF early- type galaxies from[Pasquali ctal. (2006b), 
whose SEDs were analyzed the same way as in this paper. 

main physical property tliat can be extracted from tlie 
data. We would also clarify here that we have used 'Age' 
or '<Age>' in various figures showing age distributions. 
When figures are based on single stellar population mod- 
els we use 'Age', and when figures are based on composite 
models (EXP, CSP) we use '<Age>', as they cannot give 
a single age by definition. 

4.2. Bulge Mass Estimates 

The photometry from Table [T] can be combined with 
the M/L ratios obtained from the best-fit SFH to con- 
strain the stellar mass (M^) content of the bulges. 
This M/L is derived from the composite model ob- 
tained by c ombining the simple ste llar popul ations from 
iBruzual fc Chariot (2003) using a IChabried (2003) Ini- 
tial Mass Function (IM F). I f we change the IMF from 
IChabrieil |2003 ) to Salpetol (|1955D . the stellar mass will 
increase by ~0.3 dex in log(Ms), which within the other 
errors in data and models, does not change our overall 
results. The photometry has to be corrected to take into 
account contamination from the disk, as we discuss in 
S IO We use the B/D ratio obtained from the GALFIT 
(§ I3.2p to estimate the bulge fraction of the light in the 
galaxy. The stellar mass estimates for our bulges are in 
the range of 6.5 < log(M,/MQ) < 10.0. 

Figure [Tl] shows the predicted average and RMS of 
the age distribution as a function of stellar mass and 
redshift. Over the stellar masses and redshifts probed 
in this sample we find a similar average age and scat- 
ter. This similarity could be due to two possible reasons. 
First, our sample spans a redshift range from z ~ 0.8 
to 1.3, corresponding to a difference in lookback time of 
~ 1.9 Gyr. This is comparable both to the uncertainty in 
the ag e estimate and t o the RMS of the distribution. Sec- 
ondly, Vanzell a et al.l ()2006f l found Large Scale Structure 
(LSS) in the CDF-S around z ~ 1.0 from VLT spectro- 
scopic redshifts. The redshift distribution in the HUDF 
(smaller field in the CDF-S) also show a strong peak 
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Fig. 11. — The average age [Bottom] and RMS scatter [Top] 
of the age distribution is shown with respect to stellar mass (Jcft) 
and redshift (right). Typical error bars are shown. The lookback 
time difference between redshifts z = 0.8 and 1.3 is 1.9 Gyr, i.e., 
comparable to the scatter in average ages, which explains the lack 
of a trend of average age with redshift. 

around ^~1.0. So it is possible that we may be looking 
at a smaller subset of this LSS at z~1.0. 

4.3. Disk Contamination in the GRAPES Grism 
Spectra 

The best-fit stellar population models to the GRAPES 
SEDs suggests that the late-type bulges at z ~ 1 are 
young, with an average age of ~1.3 Gyr. To better un- 
derstand this result, we first need to quantify the ef- 
fect of disk contamination in our measurements. The 
GRAPES SEDS are extracted from an aperture of rel- 
atively narrow-width aperture (5 pixels in diameter) 
around the center of each galaxy. The narrow extrac- 
tion of the grism spectra is dominated by the bulge and 
the inner disk light; since we use its 4000 A break to date 
the bulge, we need to investigate the spectral contami- 
nation due to inner disk on the estimated bulge age. We 
perform following photometric tests to understand the 
effect of the inner disk on the bulge ages. 

(1) We used the disk and bulge light profiles produced 
by GALFIT and measured their flux in a strip 5 pixel 
wide and around their common center, to estimate the 
disk and bulge light-fraction within this aperture. We 
find that the light contributed by the disk to the total 
flux in this aperture can be as high as 30%. At the same 
time, we measured the disk and bulge colors within the 
same aperture, to find that the disk and the bulge are 
similar within the photometric errors, so that the disk 
contamination in the bulge spectrum is not expected to 
dominate our estimate of the bulge age (see Appendix 
B for further discussion). This can be already be seen 
in Figure [3l where the bulge is in general 0.3-0.8 mag 
redder in {V-z') than the disk. 

(2) We compared the bulge age derived from the stellar 
population models with the color difference between two 
apertures. We measure the color difference between two 
apertures with 2.5 pixels and 5 pixels radii, equivalent 
to the narrow and wide GRAPES spectral extractions, 
respectively. The top panel of Figure [12] shows the com- 
parison between the color difference and the bulge age. 
The points on the plot are color-coded according to their 
B/D ratios, measured from GALFIT as discussed in ^ \'S.'2[ 
Blue color stands for B/D<0.5, green means 0.5<B/D<1 
and red represents B/D>1. The top panel of the Fig- 
ure [T^] does not show any major trends among age, color 
difference, and B/D ratio. Secondly, we compared the 
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Fig. 12. — Comparison between aperture colors and the best-fit 
bulge age. [Top] shows the bulge age as a function of the color 
difference between the 2.5 pixels aperture and the 5 pixels aperture. 
[Bottom] shows the bulge age as a function of the color difference 
between the 2.5 pixels aperture and the annulus defined by the 2.5 
and 5 pixels apertures. A blue colored circle stands for B/D<0.5, 
green means 0.5<B/D<1 and red represen ts B /D>1. Here B/D 
is measured in the i'-band as discussed in § 13.21 Both panels does 
not show any correlation among the age, color difference, and B/D 
ratio. 




B/D Ratio 

Fig. 13.— The Bulge-to-Disk (B/D) ratio measured with GAL- 
FIT is compared with average age [To p] a nd stellar mass [Bot- 
tom]. Similar to the trend in Figure average age docs not 
correlate with B/D either. However, the bottom panel suggests a 
correlation with stellar mass. 

bulge age to the color difference between the 2.5 pixels 
aperture and the annulus defined by the 2.5 and 5 pixels 
apertures. The bottom panel of the Figure [T^ shows this 
comparison. The points on the plot are color-coded ac- 
cording to their B/D ratios, as in the top panel. Like the 
top panel, the bottom panel of Figure [T^] does not show 
any major correlation among age, color difference, and 
B/D ratio. Finally, we directly compared the B/D ratio 
with the age and mass of the bulge. The top panel of 
Figure [T51 shows the comparison between the B/D ratios 
obtained from the GALFIT and the bulge age from the 
stellar population models for all galaxies in the sample. 
Figure [13] does not show any correlation between the age 
and the B/D ratio. Similarly, the bottom panel of Fig- 
ure [T3] shows at best a very mild correlation between the 
bulge stellar mass and the B/D ratio. 

(3) For a few sample galaxies, we extracted the spec- 
trum of their disk above and below the bulge aperture 
used to extract the bulge spectrum, at a distance of ap- 



proximately 10 pixels from the galaxy center. Similarly 
to its bulge, the inner disk also exhibits a 4000 A break 
in the spectrum whose amplitude is only slightly smaller 
(within few percents) than the bulge. This test shows 
that both bulges and inner disks are equally red/old, as 
expected from an inside-out formation scenario. We also 
fitted stellar population models using both the bulge and 
disk spectra, and analyzed how the bulge-age and metal- 
licity change as a function of disk contamination, i.e., the 
fraction of disk light added to the bulge spectrum. Our 
simulations (discussed in Appendix B) show that the de- 
termination of the bulge age is not dominated by disk 
contamination. Even when disk contamination is com- 
pletely ignored - or fully subtracted - the bulge ages do 
not change much. 

In summary, we do not detect any significant correla- 
tion between the bulge age, B/D ratio and the aperture 
color difference. We thus conclude that our estimate of 
the bulge ages is fairly robust, and that the younger age 
of the sample bulges is likely real and not due to disk 
contamination. 

5. DISCUSSION 

The ages and masses of late- type bulges are estimated 
by fitting our GRAPES SEDs with stellar population 
models. Our analysis shows that bulges in late-type 
galaxies at higher redshift (z ~ 1) appear to be rel- 
atively young, with an average age ~1.3 Gyr (6.5 < 
log(Ms/M0) < 10.0) compared to early-type galaxies at 
the same redshift. This finding appears to be indepen- 
dent of the relative amount of disk- light present, or the 
color of the underlying disk. 

Figure [14] shows the stellar masses of our bulges (filled 
and open circles) compared with the best-fit average ages 
from the CSP models discussed in § H) We also include 
the sample of e arly-type galaxies from GRAPES /HUDF 
(grey triangles: [Pasquah et all l2006bl ). whose GRAPES 
spectra were analyzed in a similar way. The early-types 
sample covers a wider range of redshifts (0.5 < z < 1.1). 
Hence, for a proper comparison, we divide both bulges 
and early-type galaxies with respect to redshift roughly 
about the median value for each subsample z ^ I for 
bulges and z ~ 0.65 for early-types. The solid (hollow) 
symbols correspond to the lower (higher) redshift bin, 
respectively. The bulges in our late-type spirals span 
a much lower range of ages, and, obviously, have lower 
st ellar masses, compared to those of early- type galaxies. 

lElmegreen et al.l (|2005l ) have classified ^900 galax- 
ies (larger than 10 pixels or 0.3") in the HUDF ac- 
cording to morphology and their photometric proper- 
ties. They find 26 9 spiral galaxies in the HUDF. Us- 
ing the lElmegreerTet al.l (j2005 ) morphological classifica- 
tions, and accurate spectro-photometric redshifts from 
I Ryan et all (|2007D . we estimate that the results in this pa- 
per represent approximately ^--^40-50% of the total late- 
type/spirals HUDF galaxy population within the magni- 
tude and redshift range used in this paper. 

Our analysis of the central and the inner disk colors 
of these galaxies (Figure [3]) and their grism spectra (as 
discussed in the Appendix B) shows that the inner disk 
and the bulge components have similar colors, and that 
the bulge ages are not significantly affected by the light 
(and stellar populations) of the underlying disk. This 
result is consistent with the idea that the inner disk of 
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Fig. 14. — Comparison between the ages of the bulges 
in this sample (black circles/lines) and early-ty pe galaxies in 
GRAP ES /HUDF (grey triangles/dashed lines; IPasQuali et al.l 
I2006bl ). The inset shows the histogram of redshifts for both sam- 
ples with ellipticals peaking around z ~ 0.65. Both samples are split 
with respect to redshift, with solid symbols representing the lower 
redshift subsample. There is a very significant difference between 
the average age of early- type galaxies and galaxy bulges. Further- 
more, the age difference is better defined for early-types, suggesting 
passive evolution for these galaxies and a more extended star for- 
mation history for the bulges. 

galax ies in general has similar c olors and age as the bulge 
(e.g.. iPeletier fc Balcellg I1996D . The effect of dust on 
these measurements should not be significant, since our 
bulge ages are based on the amplitude of the 4000 A 
observed in the GRAPES grism spectra, which is mostly 
sensitive to age and has a weak d e pende nce on dust (see 
Figure il). Also. lMacArthur etal.1 (j2004[ ) argue that dust 
is generally not a significant contributor to galaxy colors 
in low-mass/low-luminosity spiral galaxies, but is likely 
important in more massive/brighter galaxies. On the 
other hand, even if it plays an important role in this 
analysis, the inclusion of dust will make our ages even 
younger, and our result that bulges and inner-disks have 
similar dominant stellar population with an average age 
of ^^1.3 Gyr should then be viewed as an upper limit. 

We performed GALFIT decomposition on the sample 
galaxies by simultaneously fitting the bulge to a Sersic 
profile and the disk to an exponential profile. For all 
bulges in our sample, we obtained better fits using Sersic 
indices of n < l-O- Therefore, these bul ges arc disk-like 
(jKormendv fc KennicuttI 120041: 1 Athanassoulal 120051 ) and 
have radial surface brightness profi les similar to d isks. 
Similar analyse s for l ocal spirals bv lde Joii3 ()1996[ ) and 
iCourteau et al.l ()1996l ) have shown that the majority of 
bulges in late-type galaxies are better fit by exponen- 
tial profile. Our results show that a similar trend also 
exists at z ~ 1. The similarities we find in the bulge 
and the inner-disk properties (4000 A break, colors and 
profiles) could imply that these less massive, younger 
bulges at 2: 1 grow through s ecular evolution processes 
(jKormendv fc KennicuttI l2004f ). At 2; ~ 1, it is possible 
that we are seeing these galaxies still forming, and these 
"disk-like" bulges might grow from disk material or mi- 
nor mergers to become more massive bulges observed at 
present day. Disk-like pseudo-bulges can also grow by gas 
inflow an d star-for mation. Bars can drive central gas in- 
flows ( Sheth et al.|[2b05> ). and therefore, there could be a 
correlation between these disky bulges and central bars. 



iSheth et all (|2008f) find that the bar fraction m very mas- 
sive, luminous spirals is constant from 2:~0 to z~0.84, 
whereas for low-mass, blue spirals it declines significantly 
with redshift to about ~20% at z ^ 0.84, in dicatin g that 
some bars do form early enough. Elmcgre en et al.l (|2005l ) 
has morphologically classified few (~10%) of our sample 
galaxies as barred galaxies, so it will be interesting to 
investigate these late-type galaxies in future studies to 
understand this relation. 

Aperture color analysis by lEUis et al\ (|2001h for bulges 
at z<0.6 in early-type and spiral galaxies with Iab < 24 
mag found that their central colors are redder than the 
surrounding outer disk colors, but that these central col- 
ors are bluer than those of pure ellipticals at the same 
redsh ifts. As s hown in Figure [TOl our results agree 
with lEUis et al.l (|2001h . This is not perhaps surpris- 
ing, since we also select our sample based on galaxy 
total-magnitudes, with no constraints o n its bulge mag- 
nitude. In comparison, iKoo et al.l (|2005| ) select their sam- 
ple based on bulge luminosity, and they find that lumi- 
nous, high-redshift (0.73 < z < 1.04) bulges {Iab < 24 
mag) within the Groth Strip Survey are very red/old. 
They clearly show that if the bulge sample is luminous, 
then all bulges are equally red and old. In contrast, we 
show that if the bulge sample is selected without any 
constraint on the bulge magnitude, then late-type bulges 
are younger than bulges in early-type galaxies at similar 
redshifts. 

Galaxy colors and structural properties show a bi- 
modal distribution, separating into a red sequence, 
populated by early-type galaxies, and a blue "cloud" , 
populated by la te-type galaxies (|Balogh et al.l 12004 
iDriver et al.l l2006l ). Whether a galaxy resides in a red 
sequence or a bl ue cloud is also relate d to the type of 
bulge in a galaxy fPror v fc FisheH[2007f) . Figure HI shows 
this bimodal distribution. Our late-type galaxies with 
pseudo bulges lie in the bluer cloud compared to early- 
type galaxies that lie on the red color sequence. This 
shows that the processes involved in the formation of 
galactic bulges and their host galaxies are very similar. 
Observations indicate that these formation mechanisms 
depend strongly on the bulge (as well as galaxy) mass, 
and that they were active at z 0:^ 1.3. This evid e nce is 
stre ngthened by the r esults of iThomas fc DaviesI (pOOfih 
and lMacArthur et al.l (|2008( ). who find that bulges of sim- 
ilar mass have a similar evolutionary path. Possibly be- 
cause of cosmic variance, we do not detect early-type 
spirals at z ~ 1 in the HUDF. Comparing our sample of 
late-type bulges with the m assive early-typ e galaxi es at 
similar redshifts studied by iPasquali et al.l ()2006bf ). we 
confirm the existence of different evolutionary patterns 
for bulges in early- and late- type galaxies (see Figure 2]). 

Our analysis of the deepest optical survey, the HUDF, 
along with the deep unique ACS grism spectroscopy pro- 
vides the best spatial resolution at z~ 1, and yields more 
detailed insight in the process of galaxy formatio n. Mas- 
sive and luminous bulges (e.g.. IKoo et al.ir2005( ). which 
mostly reside in early-type galaxies and in earlier-type 
spiral galaxies, are old a nd formed at z>,2. The s ecular 
evolution suggested by iKormendv fc KennicuttI ()2004D 
does not play any role in their formation. Our ACS 
grism study of the HUDF here has shown that lower- 
mass bulges, which are mostly associated with later-type 
galaxies, are on average younger than one would expect 
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by letting their stellar populations passively evolve since 
their formation redshift z ~ 2. This result would point 
to secular evolution as a likely mechanism to support 
prolonged star for mation in low-mass b ulges, which is 
very different from Pasauali et al.l ()2006b) ellipticals, for 
which a quick SFH plus passive evolution can explain 
their observed SEDs. 

6. SUMMARY 

We estimated the stellar ages and masses of 34 bulges 
of late-type galaxies from the HUDF by fitting stellar 
population models to the 4000 A break observed in deep 
ACS GRAPES grism spectra. This study takes advan- 
tage of the exceptional angular resolution and depth of 
the GRAPES/HUDF data, which allow us to identify 
the small bulge components of this sample of galaxies at 
z~ 1, both on the direct and on the grism images, and to 
extract its corresponding spectrum, a method currently 
unfeasible with ground-based spectroscopy. We find that 
bulges in late- type galaxies at higher redshift (z~l) ap- 
pear to be significantly younger (with an average age of 
~1.3 Gyr) than early-type galaxies at similar redshifts. 
This finding is robust against the amount of disk-light 
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that may contaminate in part the bulge spectrum. The 
lack of a trend between average age and redshift (see 
Figure [H]) suggests star formation in bulges is extended 
over much longer times compared to early-type galaxies. 
Our results support the scenario where low-mass late- 
type bulges form through secular evolution processes. 
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APPENDIX 

EFFECT OF METALLICITY PRIORS ON BULGE AGES 

The modelling of stellar populations in this paper includes a mild prior on the distribution of bulge metallicities. For 
each choice of a star formation history (characterized by the parameters discussed in § [4]) we multiply the likelihood 
from the distribution by a Gaussian prior with mean \og{Z/ZQ) = —0.1 and RMS a{logZ/ZQ) — 0.5. The effect 
of this prior is illustrated in Figure fTSl for one of the bulges (ID 501). The marginalised distribution in metallicity, 
average age and age width (RMS) is shown for the modelling without (top) and with {bottom) the metallicity prior. 
The Gaussian prior is also shown as a dashed line in the leftmost panels. The main effect of this prior (and the 
motivation to our inclusion of this term in this paper) is that the metallicity likelihood has a monotonic increase for 
high metallicities, above solar. Hence, we decided to reduce weight for those models with very high metallicities by 
the inclusion of the prior. This decision is justified by the fact th at population synthesis m odels are not as accurate 
at super solar metallicities due to the lack of proper calibrators (jBruzual &: Charlotll2003f l. Furthermore, the range 
of metallicities for which the prior has a mild effect corresponds to the range coinmonly found in thes e systems from 
ground-based spectroscopic observations fe.g.. IShaplev et aLll2005l ; [Schiavon et al.ll2006l : iLiu et al]|200 8[). The effect of 
the prior on the age distributions is negligible (middle and rightmost panels of Figure fT5|) . For those reasons we believe 
the inclusion of the prior does not bias in a significant way our conclusions regarding stellar ages. 

But we can also show that the prior should not have a strong effect on the distribution of metallicities either. 
Figure [16] compares the predicted metallicities of our sample (histogram) . The filled circle and error bar represents 
the mean and RMS of the metallicity distribution. In comparison, the Gaussian prior (solid line) is much wider. Only 
when the prior and the predicted values have the same distribution one could suspect of strong biasing. 

Finally, in order to quantify the effect of metallicity priors on the predicted ages and metallicities for our sample, 
we show a comparison for EXP models with (black dots) and without priors (grey dots) in Figure [iTl The rightmost 
panels show histograms of the predicted ages and metallicities. As expected, the metallicities obtained without a prior 
are higher compared to models with no prior. However, taken into account reasonable uncertainties for the metallicity 
estimates extracted from unresolved spectra (of order 0.3 dex), one can say that our methodology is acceptable within 
the expected error bars. 

EFFECT OF DISK CONTAMINATION ON BULGE AGES 

In order to assess the effect of disk contamination on the age and metallicity estimates, we performed an extraction 
of the SED of the disk component for a small sub-sample of galaxies. This corresponds to the stacking of two strips 
(5 pixels wide) at an approximate distance of 10 pixels on either side of the center of the galaxy. When comparing 
the observed spectra with synthetic models, we assume that a fraction of the flux in each galaxy comes from the 
disk. Hence, for each choice of parameters, a synthetic SED is obtained, and a fraction of the disk SED is 
added to this spectrum, before performing a maximum likelihood analysis. This fraction is measured in the i'-band 
(F775W) . Figure [TH] shows the resulting ages and metallicities for four galaxies from our sample as a function of disk 
contamination. The 95% confidence levels are shown as error bars. The ages and metallicities are shown in the same 
plot, as solid and open dots, respectively. One can see that the effect of disk contamination is small, as expected from 
the weak color gradients found in the images. Th is result is consistent with the idea that the inner disk of galaxies is 
as old as the bulge fe.g.. iPeletier fc Balcellsl[l996l ). 
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log(Z/Zo) <Age>(Gyr) RMS Age (Gyp) 

Fig. 15. — One typical example (ID 501) of the likelihoods obtained for metallicity, average age and RMS of the age distribution for a 
set of EXP models with and without metallicity prior. This figure shows that without the prior we get a slightly higher metallicity but the 
average stellar age and scatter remains very similar to the modelling with a prior. 




log(Z/Ze) 

Fig. 16. — The distribution of bulge metallicities estimated from the stellar population models. The data point shows the average 
value and RMS scatter of the distribution. We overplot a mild Gaussian prior used on the metallicity with average log(Z/ZQ) = —0.1 
and RMS=0.5 dex. This prior allows for a wi de range of a verage metallicities, and is compatible with the values obtained at z ~ 1 from 
ground-based spectroscopic observations fe.g.. IShaplev et"al..,2005. : .Schiavon ct al... 2006. : .Liu et al...2008i) . 
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Fig. 17. — Comparison of bulge age distribution and bulge metallicity for EXP models with and without a metallicity prior. The age 
histogram does not show any appreciable change. The metallicity histogram without a prior (grey) peak at slightly higher metallicities — 
as illustrated in Figure [TSl — compared to the metallicity histogram (black) with prior, but it is nevertheless compatible given the typical 
uncertainty (of order 0.3 dex) in estimates of metallicity from unresolved stellar populations. 
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Fig. 18. — Effect of disk contamination on the predictions of age and metallicity. The solid (empty) dots show the estimated average age 
(metallicity) for a grid of exponentially-decaying star-formation histories. The 95% confidence levels are shown as error bars. The result is 
presented as a function of the contamination of the disk, as measured in the F775W band (see text for details). 
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TABLE 1 

Our Sample - HUDF IDs, Coordinates, Optical Magnitudes and Colors 



HUDF RA (J2000) DEC (J2000) B V i' z' {V-z')'' 



ID 








(mag) 


(mag) 


(mag) 


(mag) 


color 


501 


53. 


,1677662 


-27.8166951 


24. 


,73 


24, 


.47 


23. 


.94 


23. 


.40 


1, 


.068 


521 


53. 


,1699372 


-27.8178736 


26. 


,27 


25, 


.98 


25. 


.28 


24. 


.94 


1. 


.041 


901 


53. 


,1681788 


-27.8129432 


24. 


,44 


24. 


.03 


23. 


.35 


23. 


.06 


0. 


.970 


3048 


53. 


,1659121 


-27.7997159 


26, 


63 


26. 


.21 


25 


.40 


25. 


.17 


1 


.041 


3299 


53. 


,1923333 


-27.7978741 


25. 


,71 


25. 


,40 


25 


,00 


24. 


.59 


0, 


.817 


3372 


53. 


,1761793 


-27.7961337 


22. 


90 


22. 


,34 


21 


,64 


21. 


,25 


1 


.095 


3373 


53. 


,1668679 


-27.7976989 


24. 


,48 


24. 


,43 


23 


,94 


23. 


,77 


0, 


.662 


3613 


53. 


,1567527 


-27.7955981 


23. 


,80 


23. 


,34 


22 


,63 


22. 


,12 


1 


.219 


4072 


53. 


,1278092 


-27.7950828 


24. 


,66 


24. 


,49 


24 


,26 


23. 


,73 





.759 


4084 


53. 


,1278372 


-27.7947976 


24. 


,61 


24. 


,44 


24, 


,05 


23. 


,56 


0, 


.888 


4438 


53. 


1376781 


-27.7919373 


23. 


,58 


23. 


,13 


22 


,44 


22. 


,08 


1 


.049 


4491 


53. 


1675702 


-27.7925214 


23. 


,93 


23. 


,68 


23 


,24 


22. 


,89 





.791 


4591 


53. 


1713278 


-27.7929428 


25. 


,08 


24. 


,83 


24, 


,16 


23. 


,84 


0, 


.990 


5190 


53. 


1450905 


-27.7894219 


24. 


,22 


24. 


,00 


23 


,70 


23. 


,17 





.827 


5405 


53. 


,1606000 


-27.7897302 


26. 


,09 


25. 


,79 


25 


.12 


24. 


,64 


1 


.150 


5417 


53. 


,1661791 


-27.7875215 


23. 


,10 


22. 


,61 


21, 


.97 


21. 


,48 


1. 


.135 


5658 


53. 


,1740361 


-27.7880062 


24. 


,99 


24. 


,61 


23 


.99 


23. 


,43 


1. 


.178 


5805 


53. 


,1920649 


-27.7871824 


23. 


,91 


23. 


.58 


23. 


.04 


22. 


,57 


1. 


.011 


5989 


53. 


,1609549 


-27.7864996 


24. 


,81 


24. 


.58 


24. 


.05 


23. 


,55 


1. 


.028 


6079 


53. 


,1394080 


-27.7867760 


25. 


,76 


25. 


.35 


24. 


.78 


24. 


,12 


1. 


.222 


6785 


53. 


,1915603 


-27.7826687 


24. 


,15 


23. 


.88 


23 


.57 


22. 


,97 





.907 


6821 


53. 


1782201 


-27.7830771 


24. 


13 


23. 


,90 


23, 


,47 


23. 


,04 


0, 


.865 


7036 


53. 


1903447 


-27.7820005 


24. 


,43 


24, 


,29 


24, 


,09 


23. 


,57 


0, 


.714 


7112 


53. 


1658806 


-27.7815379 


24. 


,65 


24, 


,09 


23 


,41 


22. 


,83 


1, 


.258 


7559 


53. 


1587381 


-27.7705348 


23. 


,93 


23, 


,54 


22 


,86 


22. 


,48 


1, 


.057 


8125 


53. 


1729989 


-27.7778393 


23. 


,95 


23, 


,85 


23, 


,45 


23. 


,14 


0, 


.716 


8551 


53. 


1517635 


-27.7754183 


23. 


,90 


23, 


,42 


22, 


,71 


22, 


,25 


1, 


.173 


8585 


53. 


1479310 


-27.7739569 


22. 


,42 


22. 


,13 


21, 


,63 


21, 


,15 


0, 


.973 


9018 


53. 


1470806 


-27.7784243 


24. 


,44 


24. 


,32 


23, 


,96 


23. 


,66 





.663 


9125 


53. 


1663274 


-27.7685923 


24. 


,13 


23. 


,83 


23, 


,36 


22. 


,71 


1, 


.121 


9183 


53. 


1601763 


-27.7693039 


24. 


,86 


24. 


,80 


24, 


,40 


24. 


,18 


0, 


.611 


9341 


53. 


,1598624 


-27.7668404 


24. 


,88 


24. 


,68 


24, 


.12 


23. 


,72 


0, 


.953 


9444 


53. 


,1554162 


-27.7660748 


25. 


,81 


24. 


,91 


24, 


.04 


23. 


,34 


1, 


.572 


9759 


53. 


,1596949 


-27.7622834 


24. 


,76 


24. 


,62 


24, 


.34 


23. 


,82 





.803 



Note. — Magnitudes are in the AB photometric system, 
t Observed {V-z') color corresponds to rest-frame (U-B) color at 2~1. 
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TABLE 2 

Our Sample - IDs, Coordinates and Redshifts 



HUDF 
(ID) 


RA (J2000) 
(deg) 


DEC (J2000) 
(deg) 


phot 


VLT 


SED Fit 


501 


53.1677662 


-27.8166951 


1.07 




1.02 


521 


53.1699372 


-27.8178736 


1.04 




0.99 


901 


53.1681788 


-27.8129432 


0.94 




0.88 


3048 


53.1659121 


-27.7997159 


0.84 




0.85 


3299 


53.1923333 


-27.7978741 


1.06 


1.221 


1.20 


3372 


53.1761793 


-27.7961337 


1.04 


0.996 


0.98 


3373 


53.1668679 


-27.7976989 


0.89 




0.93 


3613 


53.1567527 


-27.7955981 


1.07 


1.097 


0.98 


4072 


53.1278092 


-27.7950828 


1.27 


1.189'' 


1.20 


4084 


53.1278372 


-27.7947976 


1.08 




1.06 


4438 


53.1376781 


-27.7919373 


1.04 


0.998 


0.97 


4491 


53.1675702 


-27.7925214 


1.05 




1.01 


4591 


53.1713278 


-27.7929428 


0.92 




0.90 


5190 


53.1450905 


-27.7894219 


1.22 


1.316 


1.18 


5405 


53.1606000 


-27.7897302 


1.07 


1.096 


0.90 


5417 


53.1661791 


-27.7875215 


1.10 


1.097 


1.03 


5658 


53.1740361 


-27.7880062 


1.07 


1.096 


1.05 


5805 


53.1920649 


-27.7871824 


1.05 




1.02 


5989 


53.1609549 


-27.7864996 


0.95 


1.135 


1.00 


6079 


53.1394080 


-27.7867760 


1.16 


1.298 


1.16 


6785 


53.1915603 


-27.7826687 


1.14 




1.29 


6821 


53.1782201 


-27.7830771 


1.07 




1.02 


7036 


53.1903447 


-27.7820005 


1.21 


0.743'* 


1.16 


7112 


53.1658806 


-27.7815379 


1.07 




1.04 


7559 


53.1587381 


-27.7705348 


0.96 




0.93 


8125 


53.1729989 


-27.7778393 


1.05 




1.04 


8551 


53.1517635 


-27.7754183 


1.05^= 


1.047 


1.05 


8585 


53.1479310 


-27.7739569 


0.97 


1.088 


1.05 


9018 


53.1470806 


-27.7784243 


0.95 




0.99 


9125 


53.1663274 


-27.7685923 


1.20 


1.295 


1.21 


9183 


53.1601763 


-27.7693039 


0.94 




0.94 


9341 


53.1598624 


-27.7668404 


1.03 




1.01 


9444 


53.1554162 


-27.7660748 


1.07 


1.096 


1.02 


9759 


53.1596949 


-27.7622834 


1.34 




1.22 



^ From GRAPES photometric redsiiift catalog lIRvan et al.ll2007D . 

^ From GOODS-IVIUSIC catalog lIGrazian et al.ll2006D and lVanzella et all 

| [200l) . 

'^ From GRAPES SED fitting (§ g}, which is dominated by the 4000 A 
break fitting. 

Quality flag on these spectroscopic redshifts is poor. 



